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Abstract 

We describe the general setting for the optical Aharonov-Bohm 
effect based on the inverse problem of the identification of the coef- 
ficients of the governing hyperbolic equation by the boundary meas- 
rements. We interpret the inverse problem result as a possibility in 
principle to detect the optical Aharonov-Bohm effect by the boundary 
measurements. 

1 Introduction. 

In this section we will review the quantum mechanical Aharonov-Bohm (AB) 
effect (c.f. [AB], [WY], [OP], [W], [E4]). 

Let r2 be a smooth bounded domain in R" having the form Q = Qq\ 
lJY=i^j, where Qq is a simply-connected domain and Qj,l < j < m, are 
smooth domain called obstacles. We assume that Qj C for 1 < j < m, 
and Qj fl fi^ = when j k, 1 < j,k < m. 

Consider the stationary Schrodinger equation in Q with magnetic poten- 
tial A{x) = (v4i(x), ...,An{x)) and electric potential V{x): 

(1.1) ^i(^)) u{x) + V{x)u{x) = k^u{x), 

describing the nonrelativistic quantum electron in the classical electromag- 
netic field. We assume that 

(1.2) u\dn, =0, 1 < J < m, 
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i.e. Qj are unpenetrable for the electron, and 



(1-3) u\ong = f{x). 

Let A.{k)f be the Dirichlet-to-Neumann (DN) operator on 9^0, i-e- 



:i.4) A(^k)f=^^-^iA.Jy)u 



where u{x) is the solution of (11. II) . (II. 2p . (I1.3P and ly is the unit outward 
normal vector at x G dQo- 

Denote by G{Q) the group of all complex- valued C°°{Q) functions c{x) 
in Q such that |c(x)| = 1. 

If c{x) G G{Q) and u' = c~^{x)u{x) then u' satisfies the Schrodinger 
equation of the forn (II. ip with A{x), V{x) replaced by A'{x), V'{x), where 

(1.5) A'.{x) = Aj{x)~ic-\x) — , l<j<n, 

V'{x) = V{x). 

We shall call the electromagnetic potentials A'{x),V'{x) and A{x),V{x) 
gauge equivalent. We also call the DN operators A{k) and A'{k), corre- 
sponding to A{x),V{x) and A'{x),V'{x), respectively, gauge equivalent if 
there exists c(x) G G{^1) such that 

A\k) = c^'Aik)co, 

where Cq is the restriction of c{x) to dflo. 

Let B{x) = curl A{x) or, equivalently, B = dA, where A = J2]=i Aj{x)dxj, 
be the magnetic field in Q. It follows from (ll.Sp that 

B{x) = B\x) in n 



if A{x) and A'{x) are gauge equivalent. If Vt is simply-connected then the 
inverse is true: B{x) = B'{x) in Q implies that A{x) and A'{x) are gauge 
equivalent. When Q is not simply-connected this is not true anymore. It was 
shown in the seminal paper of Aharonov and Bohm [AB] that if curl A = 
curl A' = 0, but A'{x) and A{x) belong to distinct gauge equivalent classes, 
they have a different physical impact that is detectable in the experiments. 
This fact is called the Aharonov-Bohm effect. 
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An important description of gauge equivalence classes was given by Wu 
and Yang [WY] : 

Let 7 be any closed path in Q. It is easy to see that A{x) and A'{x) 
belong to the same gauge equivalent class iff 



where p e Z. 

In the original paper [AB] Aharonov and Bahm consider the case of one 
obstacle Qi in and the magnetic field confined to Qi. Then J^A-dx = a 
is the magnetic flux and a is independent of any simple path 7 encircling Qi. 
The quantity e*" that determines the gauge equivalence class of A{x) was 
measured in this experiment. If a 7^ 27rp,p G Z, then the gauge equivalence 
class of A{x) is nonzero despite the fact that i? = in = \ ^i- 

Consider now the case of several obstacles fli, ...,Qrn- Suppose that the 
magnetic field is hidden inside each of these obstacles. Let ak = Xy^. ^ " dx 
be the magnetic fluxes, where 7^ encircles flk only. Suppose that some of 
1^ are not integers and YlT=i = 0, i.e. the total magnetic flux is zero. 
In this case the gauge equivalence classes are determined by m parameters 
giQfc ^1 < k < m, however the AB experiment will not find a gauge equivalent 
class different from zero. To identify an arbitrary gauge equivalence class one 
needs to use broken rays (i.e. the rays reflected at the obstacles) belonging 
to the base of the homotopy group of fl (c.f. [E5], page 1512). 

It is necessary to perform at least m AB type experiments to determine 
all e*"*, 1 < k < m. When B{x) = curl A is not zero in Q it is not enough 
to perform a finite number of AB type experiments to identify the gauge 
equivalence class of A. Therefore the following question arises: Is it possible 
by the measurements on the boundary dflo to detect the difference in the 
gauge equivalence classes of A{x) and A'{x)7 The answer to this question is 
affirmative, and it is given by the following theorem (c.f. [E4], [W], [N], [KL] 
and further references there): 

Theorem 1.1. Consider two boundary value problems U.l\} . U.3\} cor- 

responding to electromagnetic potentials A{x),V{x) and A'{x),V'{x). Then 



(1.6) 




for all paths 7 in fi, or, equivalently. 



(1.7) 
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A(x), V{x) and A'{x), V'{x) belong to the same gauge equivalence class iff the 
DN operators A{k) and A'{k) are gauge equivalent for all k. 

We consider each boundary measurement as an experiment. The Theorem 
11.11 asserts that the boundary measurements are able to identify an arbitrary 
gauge equivalence class. We interpret this theorem as a confirmation of the 
Aharonov-Bohm effect. □ 

In §2 we develop the same approach in the case of the optical Aharonov- 
Bohm effect, and we shall formulate the unique identification theorem for 
the optical AB effect. In §3 we prove the main unique identification the- 
orem (Theorem 12.31) . Our approach to the hyperbolic inverse problems is 
based on a modification of the BC-method given in [El], [E2]. The power- 
ful BC-method was discovered by M.Belishev and extended by M.Belishev, 
Y.Kurylev, M.Lassas and others (c.f. [B], [KKL], [KL] and additional refer- 
ences there). An important part of the BC-method is the unique continuation 
theorem by Tataru[T]. The approach of [El], [E2] allows one to consider new 
problems that were not accessible by the BC-method as the inverse hyperbolic 
problems with time dependent coefficients (see [E3]). The inverse problem 
results of this paper are also new. 

2 The optical Aharonov-Bohm effect. 

In this section we consider hyperbolic (wave) equation of the form: 

where x = (xi, x„) G Q, xq is the time variable, ([5f-'^]"^^o)~"^ is the pseudo- 
Riemannian metric tensor with Minkowsky signature, i.e. the quadratic form 
Y.lk=o9^''(^)^j^k has the signature (1, -1, -1), g{x) = {dety''])-\ We 
assume that g^^{x) are smooth in Q and independent of Xq. 
We make two additional assumptions: 

(2.2) (1,0,..., 0) is a time-like direction, i.e. g^^{x) > 0, x G fi. 
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and 
(2.3) 

n 

The plane = intersects the cone g^^{x)^j^k = at 

i,fc=0 

(ei,e2,...,en) = (0,0,...,0) only, 

n 

i.e. the form — g^^{x)^'j^k is positive definite, x eVL. 

j,k=i 

The important physical example of equation of form fl2.1l) is the equation of 
the propagation of light in the moving medium. Here the tensor g^^{x) has 
the following form (see Gordon (1923), [NVV], [LPl]): 

(2.4) g^^ = + {n\x) - l^u^ < j, k < n, n = 3, 

when [r/-'^]""^ is the Lorentz metric tensor, rj^'^ = 0, when j ^ k, r]^^ = 
1, T]^^ = —1 for 1 < i < Xq = ct, n{x) = ^Je{x)^{x) is the refraction 
index, {yP , u^, v? , v?) is the four-velocity of the medium fiow, (w", v} , u^, u^) = 
(1 — ■^)^2(1, w{x) = {wi,W2,Ws) is the velocity of the flow (c.f. [LP], 
[LP1],'[LP2]). ' 

In the case of slowly moving medium one drops the terms of order (M)2 
(c.f. [LPl], [LP2], [CFM]). Then the metric of the slowly moving medium 
has the form: 

(2.5) g^^ = rf'' for l<j,k<n, 
g'' = n\x), g'^ = g^' = v,{x)'=' (n'-l)^, l<j<n,n = 3, 

and the corresponding equation is 

We shall also consider in addition to the equation (12.61) the following 
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equation: 
(2.7) 

2, .S^U 1 f d , \ 9 \ / , / M / / \ \ 

where Vj{x) and v?{x) are the same as in (12. 5p . 1 < i < n. 

Equation (12. 7p differs from the equation (12.61) by the term ^"=1 t'|(x)|^. 

Since f| = 0((^)^) the equation (12. 7p also describes the propagation of hght 
in the slowly moving medium. We consider (12. 7p to have a closer analogy 
with the quantum mechanical AB effect, although the addition of extra terms 
affects the uniqueness of the inverse problem (compare Theorems 12.11 and 
12.20 . Note that the nonuniqueness is of the first order in — (see Theorem 

EI]). 

We consider the initial-boundary value problem for (12.61) and (12. 7p in the 
infinite cylinder Vt x (—00, +00), where VL is the same domain as in §1: 

(2.8) u{xq^x) = Q for <^ 0, 



(2.9) M(a^o,a;)|aQ.x(-oo,+oo) = 0, 1 < i < m, 

M(a;o,a^)|9Cox(-oo,+oo) = f{xo,x'), x' E dQo, 

where f{xo, x') has a compact support on OVLq x (—00, +00) 
Denote by A the hyperbolic DN operator: 



, du , , du 
(2.10) Kf=[ [VP 



1 9f2o X (—00, +00) ) 



dp dxQ ^ 

where, as in §1, is the external unit normal to dQo- 

In studying the equation (12.70 we shall use the following change of vari- 
ables in n X (—00, +00): 

(2.11) xq = xq + a{x), Xj = Xj, I < j < n, 

where a{x) G C°°(f2), a{x) = on dQo. If u{xo,x) is u{xo,x) in new coordi- 
nates, then u{xo, x) also satisfies an equation of the form (12. 7p : 

(2.12) Mt/f.2.,.^Mi:o,x) 



dxl 



— ; ( t: Vjix)^^] \AM^\ ( t: Vj{x)^z- ] "U = 0, 

j^i VWM V^^J- dxj ^ \dx^ 'dxoj 
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where fi{x) = n{x), Vj{x) is replaced by 

(2.13) Vj{x) = Vj{x) - a^^{x), 1 <j <n. 

We assume that < n^{x) to preserve the hyperbohcity of fl2.12p . 

Note that 

(2.14) M = for xo < 
and 

(2.15) u\dn,x(-oo,+oo) = 0, I < j <m, 

"^liSCoX (— oo,+oo) 

= f{xo, x'), 

where /(xq, 2;') = /(xq, x') since a = on OQq 

We shall say that Vj, I < j < n, and Vj, I < j < n, belong to the same 
equivalence class if (12.131) holds. 

If v{x) = {vi{x), ...,Vn{x)) and v{x) = {vi,...,Vn) belong to the same 
equivalence class then 

(2.16) / V ■ dx — / V ■ dx = 

J 'y J 'y 

for all closed paths in Q since J2]=i Xy '^x^dxj = 0. 

It is easy to see that if v and v belong to the same equivalence class then 
the DN operators A and A are equal on OQq x (—00, +00). A nontrivial fact 
is that the inverse statement is also true. The following unique identification 
theorem holds: 



Theorem 2.1. Let Lu = 0, Lu = be equations of the form ^2. 7| ), \2.12^) 



in domains Vt,Vt = ^Iq \ U^^i^j, respectively. Suppose that the DN op- 
erators A and A are equal on OQq x (—00, +00) for all f G C^{dflo x 
(— cx),+oo)). Then the ^2 = ^2, h{x) = n{x) and the corresponding velocity 
flows v{x), v{x) belong to the same equivalent class, i.e. Ii2.13\) holds for some 
aix)EC^m, a(x)|ano = 0. 

Note that we did not assume apriori that the number of obstacles rh in 
Cl and their location are the same as in Q. 

A consequence of Theorem 12.11 is that boundary measurements on dflo x 
(— cxD, +CX3) uniquely determine the integrals J^v ■ dx for all paths 7 in f2. 
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As in §1 we view the optical Aharonov-Bohm effect as the fact that the dif- 
ferent equivalence classes of the velocity flow have different physical impacts. 
Theorem 12.11 confirms that the boundary measurement (experiments) allow 
one to distinguish different equivalence classes, i.e. to detect the Aharonov- 
Bohm effect. 

Remark 2.1 There is a difference between the optical Aharonov-Bohm 
effect and the quantum mechanical AB effect. In the case of the optical AB 
effect the boundary measurements allow one to recover J^v-dx. In the case of 
the quantum mechanical AB effect we can recover only J A- dx (mod 27rp), 
peZ. ^ □ 

Let u{k, x)) = )e (ixo be the Fourier-Laplace transform of 

u{xo,x) in xq, or let e^^^^^u{k,x) be a monochromatic wave. Then u{k,x) 
satisfies the Schrodinger equation: 

(2.17) - en\x)u(k, x) - V , ( ikvJx] 

d_ 



^^^)\\^~ (a;) \u{k,x)=^ 



with the boundary conditions 

(2.18) {i(A:,x)|an, =0, 

u{k,x)\dQ.^ = f{k,x'). 

Now kv{x) plays the role of the vector potential and it depends on k. Note 
also that the Fourier-Laplace image T of the transformation (12. lip is the 
multiplication by e^^°-^^) ^ i.e. T is a gauge transformation depending on pa- 
rameter k. 

When VL is multi-connected one can expect that the Aharonov-Bohm ef- 
fect takes place for fl2.17p . This problem was studied in optics (c.f. [LP], 
[LPl], [LP2], [CFM]). An analogous problem was considered for the water 
waves and for the acoustic waves in [BCLUW], [RdeRTF], [VMCL]. 

These authors considered the case of one obstacle Vti C and irrota- 
tional flow in f2o \ ^^i. Performing an Aharonov-Bohm type experiment they 
measured exp(z J^v ■ dx) as in the quantum mechanical AB effect. Since 
such experiments are based on the geometric optics considerations it was as- 
sumed that the light rays are straight lines and kv{x) is not large. A natural 
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question arises whether some form of the AB effect takes place when these 
conditions are not satisfied. 

Note that a rigorous geometric optics approach when k —>■ oo for the 
equation fl2.17p is more dehcate than for the equation (11. ip . In particular, 
the eiconal equation depends on v{x) and the light rays are not the straight 
lines. 

Remark 2.2 Let curl v = in = {VIq \ i^i)) C R^. In this case the 
equivalence class of v{x) depends only on one parameter a = J^v ■ dx, where 
7 encircles Qi. There is a simple solution of the inverse problem in this case 
that does not use neither the geometric optics nor the Theorem 12. H 

Let v{x) and v{x) be two irrotational velocity flows in Qq \ fii. Consider 
two Schrodinger equations of the form f l2.17p in = \ fii assuming that 
77,(2;) = n{x) in Q and A(fc) = A(A;) on dQo for some fixed k. It was shown 
in [NSU], [ER], using the parametrix of the DN operators, that v ■ t = v ■ r 
on dVto, where r(x) is the tangent vector to dVLo at a; e SVLq. It follows from 
V ■ T = V ■ T on OVLq that a = J^^^ v ■ dx = Jg^^ v ■ dx. Since v and v are 
irrotational this implies that there exists a{x) E C°°(fi) such that v — v = 
Since |^ ■ r = on dilo we get that a|aQg = = const. Replacing a{x) by 
a{x) — Oo we obtain that v and v belong to the same equivalence class. 

Similar arguments apply in the case of equations (12. 6p and (12. ip with 
the metric (12. 4p . Using the parametrix of the DN operator we can recover 
the restriction of the metric to OQq (c.f. [LU] or [El], Remark 2.2). In 
particular, we can determine w{x) ■ t{x) on dQo. Therefore we can recover 
a = Jg^^ w{x) ■ dx. In the case of irrotational flow and one obstacle a is the 

same for any simple path in Q = fio \ ^i- D 
We shall investigate now the inverse problem for the equation (12.60 . The 
case of the equation (12.10 with the metric (12.40 will be studied in another 
paper. 

Theorem 2.2. Consider two initial- boundary value problems in domains Qx 
{—00, +00) and Cl X (—00, +00) for operators of the form h2. corresponding 
to the metric tensors [g^^{x)Y'^^ [g^^{xy\~^ of the form / (^.5|) . respectively. 
Assume that the DN operators A and K, corresponding to L and L are equal 
on OVLq X (—00, +00). Assume also that there exists an open dense set O gVL 
such that the velocity flow v{x) = {vi, ...,£'„) does not vanish on O. Then 

Q = Q, h{x) = n{x), v{x) = v{x), 1 < J < 

unless v{x) is a gradient flow, i.e. there exists b{x) G C°°(fi) such that 
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v{x) = ^ and b(x) = on O^Iq. In the case of the gradient flow there are 
two solutions v{x) = v{x) and v{x) = —v{x). 

The proofs of Theorems 12.11 and 12.21 will be given in the end of this 
section. □ 

Now we shall consider the general case of the initial-boundary value prob- 
lem (I23D, (E9]) for the equation (O). The DN operator for (El]) has the 
following form: 



(2.19) A/= J29'\^) — ''k 



du 

' dxj 

j,k=0 



n 



p,r=l 



1 9f2o X oo,+oo) ) 



where v is the unit normal as in (12.101) . 

Consider the diffeomorphism of the form: 

(2.20) xo = + a(x), 

X = ip{x), 



where a(x)|an„ = and f{x) is a diffeomorphism of Q onto Q, where f2 is a 

domain of the form = ^lo\ WJ^i^lj and ip = I on O^Iq. Note that (12.201) 
transforms (12. ip into the equation of the same form. More precisely, (12.11) 
has the following form in {xo,x) coordinates: 



where 



(2.22) W'^ix)] = Jix)y\x)]J^ix), 
g{x) = {det[g^\x)])-\ 0<j,k<n, 

J{x) is the Jacobi matrix of (I2.20p . 

Theorem 2.3. Consider equations ^2. 1\) and ^2.21\) in domains fix (— oo, +oo) 
and fix (—oo, +oo), respectively, with initial-boundary conditions ^2.8i) . ^2. 9\) 
and l\2.14\ l> i^2.15\) . respectively, where f = f- Assume that the DN operators 
A and A are equal on dQo x {—oo, +oo) and the conditions Ii2.2\) . h2.3\) hold 
for L and L. Then there exists a map i/j of the form l{2.2(]\) such that 

(2.23) ipoL = L m fix(-oo,+oo). 
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Note that fl2.23l) is equivalent to fl2.22l) . Note also that since is a 
diffeomorphism of Q onto Cl, we have that rh = m and dflj are diffeomorphic 
to dClj, 1 < j < m. 

The proof of Theorem 12.31 will be given in §3. 

Remark 2.3. Making the Fourier-Laplace transform in (12. ip we obtain 

d 

(2.24) L(ik,—)u(k,x) = in Q, 

ox 

where L{-^, ^) is the operator (12. ip . Let A(A;) be the Fourier-Laplace im- 
age of the DN operator (I2.19p . Using well known estimates for the initial- 
boundary value problem (12. ip . (12. Sp . (12. 9p one can prove that the hyperbolic 
DN operator (I2.19P on OQq x (— oo, +oo) uniquely determines the DN opera- 
tor A(A;) for the elliptic boundary value problem (I2.24p . (I2.18P and vice versa 
(see, for example, [KKLM]). Here k E C \ Z, where Z is a discrete set. 

Suppose g^'^ — r]^^ = 0, when |x| > R, and suppose that fio D {a: : 
\x\ < R}. It is well known that A{k) given on dQo for fixed k = ko 
uniquely determines the scattering amplitude a{9,uj,k) for k = k^ and any 
9 G S"'~^, ijj G iS""^, and vice versa (see, for example, the recent work [OD] 
and additional references there). 

Therefore one can consider the inverse scattering problem for (I2.24p in 
R"' instead of the inverse boundary value problem for (I2.24p . (I2.18p . In 
the case when there is no obstacles and the principal part of (I2.24p is the 
Laplacian, such inverse problems were studied for r?, > 3 and fixed k (see, 
for example, [NSU] and [ERl], where the case of exponentially decreasing 
electromagnetic potentials was considered). When obstacles are present or 
when the metric is not Euclidean the hyperbolic inverse problem approach 
is much more powerful. □ 

We shall show now how Theorem 12.31 implies Theorem 12.11 and Theorem 



Proof of Theorem 2.1 Consider two equations of the form (12. 7p and 
f lZT^ . i.e. g^^ = -6jk in (EU) and g^'' = -6jk in dMH), I < j,k <n. We 
assume that A = A on dQo x (— oo, +oo). It follows from Theorem 12.31 that 
there exists a map ip of the form (I2.20p such that (I2.22p holds. It follows 
from fl2:22|) that 

(2.25) 9''=y9''^^, l<J,k<n, 



p,r=l 
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where 



(2.26) v9 = (v9i, v9„) = / ondQo- 

Since g^'' = —6jk, g^^ = —Spr we have that ^j{x) = Xj, I < j < n, is 
the solution of fl2.25p and the uniqueness of the Cauchy problem fl2.25p . 
f l2.26p implies that ip = I is the only solution of fl2.25p . (12.261) . Therefore the 
map (I2.20p reduces to the map (12.111) . This implies that Cl = Q. Making the 
change of variables (12.111) with the same a{x) as in (12.201) we get two identical 
operators. Therefore (12.131) holds and, subseqently, n{x) = n{x). □ 
Proof of Theorem 2.2. It follows from Theorem 12.31 that there exists 
a map of the form fl2:20|) such that (12:221) holds. Let [gjkix)] = [g^''{x)]-^, 
[gjkix]] = [g^^{x)]-\ Then fl^:^ is equivalent to 

n n 

(2.27) ^ gjk{x)dxjdxk = ^ gjk{x)dxjdxk, 

j,k=0 j,k=0 

where {xo,x) are related to {xo,x) by (I2.20p . Note that (c.f. [LPl]) 

(2.28) ^00 = n~~'^{x), gjk = Sjk for l<j,k<n, 
9oj = 9jo = -{n {x)-l) = n {x}Vj{x), 



and gjk have a similar form. Here Vj{x) is the same as in (12. 5p . Since 
gj^ = gi^ = —6jk for 1 < j, k < n, we have, as in the proof of Theorem 12. 1[ 
that X = (p{x) = X. Therefore Vt = VL. Note that 

(2.29) dxo = dxo + S^^^^^dx^. 

^ OXj 

Substitute (I2.29P into (I2.27p . Taking into account that n, 
and that dxo, dxi, dxn are arbitrary, we get from (12.270 and (I2.29P : 

(2.30) h-^{x) = n-^{x), 

(2.31) 2n-'^{x)a^^ + 2h-^{x)vj{x) = 2n''^{x)vj{x), I < j < n, 



(2.32) h '^(x)C^^ax^dxjY + 2(^^n ^(x)vj{x)dxj)C^^axjdxj) = 0. 
j=i j=i j=i 
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It follows from fl2.3Up that n{x) = n{x). Multiplying (12.31 1) by n^{x) we get 



(2.33) Vj{x) + aj:^{x) = Vj{x), 1 < j < n. 

If there exists x such that not all ttx^ix) =0, 1 < j < n, then we can cancel 
"'~^(^) X]j=i '^xj {x)dxj in fl2.32l) and get X]j=i {x)dxj + 2 X]j=i Vj{x)dxj = 
0, i.e. 

(2.34) a^^.(x) + 2vj{x) = 0, 1 < j < n, 

since dxj, 1 < j < n, are arbitrary. Comparing (12.331) and (12.341) we get 

y.(^x) = -Vj{x), l<j<n, 

when 1^ = (^, ^^) 7^ 0. In the case when v{x) 7^ on an open dense 
set in Q we have, by connectness of Q and by the continuity, that (12.341) 
holds on Q if a{x) ^ 0. Therefore v{x) is a gradient flow, v{x) = —v{x) 
is also a solution of the inverse problem, and it is the only solution except 
the trivial solution v{x) = v{x) that corresponds to a{x) = 0. Note that 
the boundary measurements can not distinguish between these two solutions 
v{x) and —v{x). If -0 = on an open set in then (I2.32p implies that |f = 
on this set. In such case there can be more than two solutions of the inverse 
problem. For example, if v{x) is a gradient flow, i.e -0 = b{x) = on d^lo 
and the closure of the set {x E fl : b{x) 7^ 0} is not connected, then there 
exists at least four solutions of the inverse problem. 

If v{x) and v{x) are any two solutions of the inverse problem then (I2.33p 
implies that j[^v{x)-dx = j]^v{x) ■ dx for any '-f in Q. Therefore the boundary 
measurements uniquely determine J^v{x) ■ dx. This fact can be considered 
as an analogue of the Aharonov-Bohm effect. 

3 The proof of the main theorem. 

As in [El] we start the proof of Theorem 12.31 with the introduction of a 
convenient system of cooordinates that simplifies the equation. 

Let Uq be a neighborhood of some part T of O^Iq and let ( X , x^) be a 
system of coordinates in Uq such that x„ = is the equation of dflo H Uq. 
Let T be small. 
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Denote by ip"^, the solutions of the eiconal equations in Uq 

n 

(3.1) E 9''{x)^|;iixo,x)^|^^^{xo,x) = 0, 

j,k=0 

such that 

(3.2) ip'^ = Xq when Xn = 0, 

— T — xq when Xn = 0, 



Solutions '0^(xo,a;) exist for < x„ < 5 where (5 is small. We assume that 
surfaces ij}^ = and -0" = intersect when a;„ < 6. 

In the case when g^^{x) are independent of a;o we have 

(3.4) ^+ = xo + ¥?+(x), 

= r - xo + ^~{x), 



where ^ {x) satisfy 



g'\x) ± 2 J] + g^\x)^%V% = 0, 

j=l j,k=l 



±1 _n ±1 _ T^°"(3:) + V(g°"(x))^ - ffQO(a;"y^ , 
Denote by </?p(x) the solutions of 

n 

(3.5) J]^^'=(x)V'-.(^p.,(x)=0 
with the initial conditions 

(3.6) </?pU„=o = 1 < p < n - 1. 
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Note that ^Ppxo = 0, ip^^^ = — 1. Therefore we have 



XI Vx, V'pxfc - XI Vpx, =0, 1 < p < n - 1. 

j,k=l j=l 

Make the following change of variables in Uq x [0,T]: 

(3.7) s = 'iIj~^{xo,x) = xq + (p~^{x), 

T = i/j'ixo, x) = T - xq + ^^(x), 
yj = ipj{x), l<j<n-l. 

We shall call (s, r, y') the Goursat coordinates. 

Let u{s,T,y') = u{xo,x). Then u{s,T,y') satisfies the equation 



Lu 



(3.8) 

de_f 2 d 

1^ 



/ — 

9^'~{s,r, y')^\g\ — 



n-l 



2 d_ 



-du\ 



n-l 



2 a 
2 d 



9 [s,r,y )^/\g\ — 



^+'^(s,r, y ) VW\^ 



n-l 

+ E 



1 a 



9''\s, T, y' 



j- (9n 

dyk 



The terms containing vanished because of (13. ip . (13. 5p . Here 

(3.9) ^=(-4(^+-)-^det[^^-1^,^i,)"'- 
It follows from fl3.7p that 

s + T - T = (p'^{x) + (/7~(a;), 

s — T + T = 2xo + ^p^{x) — ^p~{x). 
Denote (c.f. [El], (2.23) ) 

(3-10) yn = ^ = ^ , 



0. 



S-T + T 



Xq + 



^ - ip {X) 



2 2 

Vj = Vjix), 1 < j < n - 1. 
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We shall also use the coordinates fl3.10l) . 

Note that (^'^ = Lp~ = when Xn = 0- Therefore the map fl3.1UI) is the 
identity on x„ = 0. 

Since Ug = |('Ujyo ^ "^r = ~|(%o + Uy„) the equation fl3.8p has the 

following form in {yo,y',yn) coordinates 



(3.11) 




We used above that g^^, g^, g^'^ depend on (?/', yn) and do not depend on yo. 
Divide (13111) by g^^~ . 
As in [El] put 

(3.12) w'= |^|3(^+'-)k. 

Then u' will be the solution of the equation 

("3) (,S'^) 




where gl^ = {g+' ) ^g^'', g^^ = -g^^ = {g+' ) ^g^'^, I < j,k < n - 1, Vi 
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has a form similar to (2.8) in [El]: 



(3.14) 




, dA dA 



9o 



dA dA 
dyj dyk 



jr[\dyn\ dyjj dyj \ dynj dyjdynj' 
where A = ln[(^+'")5|^|l] = Hj^gf), 9i = {detW%;,l,)-\ = -1 (c.f. 



(ESD and (1312D ). 

Note that Liis formally self- adjoint. The DN operator Ai corresponding 
to Li has the following form: 



where /' = u'\y^=Q. It follows from the Remark 2.2 in [El] that = 
11 - 

{g^'~)^\g\^ = -^di and its derivatives on y„ = can be determined by 
the DN operator A of L. Therefore the DN operator Ai of Li is determined 
by the DN operator of L (c.f. [El], (2.9)-(2.12)). 

Introduce notations similar to [El], p. 819. Let V C F^^^ C F^^^ C 
Uq n OVLq. Denote by Dj^Q-,^ < J < 2,0 < Sq < T, the forward domain 

— (i) - 
of influence of F x [so,T] in the half-space > 0. Denote by the 

backward domain of influence of F x [0, T] for y„ > 0. 

Let YjsQ,so G [0,T),1 < j < 2, be the intersection of DjsQ with the 
plane T — yn — yo = 0. Denote by Xjs^ the part of Djsq below Yjsq. Let 
ZjsQ = dXjsQ \ (Yjso U {yn = 0}). We assume that X20 H C Uq and X2Q 
does not intersect dQ for > 0. We shall call Djso fl Dj the double cone 

of influence of F^-'^ x [scT]. Denote by Rjso the intersection of Dj^g fl Dj 
with YjsQ. 

We shall assume that F^^), 1 < j < 2, are such that Dio fl dQa C F^^) x 



Let Qj be the rectangle in the plane r = : Qj = {{s,T,y') : r = 

0,0 < s < T,y' E T )}. Note that Qj is the intersection of Dj with the 
plane r = 0. Therefore Rjsg is the intersection of YjsQ with Qj,j = 1, 2. Note 



(3.15) 




[o,r]. 
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also that if {s,0,y') G Yj^^ then the hne segment {s,0,y'),s < s < t, also 
belongs to Yjsg. Later we shall introduce one more set C r^^^ X [o,r] 
and assume that Fx [0,T] C Riq. We shall refer to this assumption and to 
the assumptions in the preceeding paragraphs as the geometric assumptions. 
These assumptions can be always satisfied if T is small. 

The following theorem is a generalization of Lemma 2.1 in [El]: 

Theorem 3.1. Let L^^^ and L^'^^ be two operators of the form Ii3.11\) and 

let A*^*) be the corresponding DN operators. Assume that A^^^ = A*^^-* on 
r*^^) X (0, T) and that the geometric assumptions are satisfied. Then there 
exists changes of variables yo = yo,yn = ymV' = ct^^Kum u') i i = li2, such 
that a'-^\0,y') = a'-^KO,y') = y' and = L^^) when y' G F, y„ e [0,f]. 
Here L^^'^ are differential operators L^*^ in the coordinates {yQ^iimy')- 

Many parts of the proof of Theorem 13.11 are the same as in Lemma 2.1 in 
[El]. We shall skip the proofs in such cases and concentrate only on the new 
elements. 

We shall start with the derivation of Green's formulas analogous to for- 
mulas (2.33) and (2.24) in [El]. 

Consider the following initial-boundary value problem for Li. 

LiU = for yn > 0, 

u = Uy^ = for yo = 0, y.^ > 0, 

where supp / C F^^^ x (0, T], F^^) cUoH {y^ = 0}. 
Let V be such that 

Llv = 0, yn> 0, 
V = Vy^ = when yo = 0, y„ > 0, 

v\y„=o = g, supp ^ c F X (o,r]. 

We have 

= {Liu, v) — {u, Llv), 
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where (m, v) = J^^^ u{yo, y)v{yQ, y)dyody'dyn. Integrating by parts we get 



Tl— 1 r\ r\ Tl — 1 r\ r\ 

O no OU „ O nn OU 



y„=o dyk A„=o^ ^l/fc 



We used here that u,v vanish on Z20. Note that other terms in Li are the 
same as in [El], formula (2.33). Therefore integrating these terms by parts 
as in [El], (2.33), and combining with fl3.16p we get the folllowing Green's 
formula: 

(3.17) / (% - u^) dsdy' = - [ {A^fg - fA^)dy'dyo, 

where Ai is the DN operator f l3.15p . Note that = Li in our case. Therefore 
the left hand side of (13.171) is determined by the boundary data. 

Now we shall derive another Green's formula similar to (2.24) in [El]. 
Consider 

dv ^ , du 

3.18 o = ^L^u,—) + i—,L,v). 

oyo oyo 

Integrating by parts in yj and s we get 



n— 1 Q n— 1 



j=l j=l -^J 



+ / J^Soj^'^dy'dyo. 

Jy„=o oyj 
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Now integrate by parts in ?/o and then again in s and yj. We get 

(3.19) (2E4«s^|.«») + pE|9?|^.«j 

j=i j=i 

Tl 1 r\ r\ 1 r-i r\ 



The remaining terms in (13.181) are the same as in [El], formulas (2.18)-(2.25). 

Therefore, combining all terms after the integration by parts we get (c.f. 
[El], (2.25)): 

(3.20) = (LiM, + {uy,,,Lxv) 

= Q{u,v)+Ao{f,g), 



where 



Y20 



(3.21) Q{u,v) = 

4„ - _ V a^^—'?^ -2 V fa'^—'^+a'^ ——) + V.uv 



dy'ds 



and 



(3.22) Ao(/, g) = / (Ai/^ + fy,Aig)dy'dyo. 

Jr(2)x[o,T] 

Again in the derivation of (13.201) we used that u = v = on Z20. 

We shall show now that the "ellipticity" condition (12. 3p . i.e. that the 
reduced quadratic form is negative definite, implies that Q{u, v) is positive 
definite. Note that the map of the form (13.71) and, consequently, the map 
(I3.10p . preserves the ellipticity condition. 

The reduced quadratic form in (I3.13P has the form: 

n— 1 n— 1 

(3.23) Yl 9l\x)i,i, -en-2j2 9o^^^n. 

j,k=l j=l 
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The "elhpticity" condition fl2.3p imphes that fl3.23p is negative definite. Re- 
placing in the complexification of fl3.23l) by 2us and by —Uy. , 1 < j < 
n — 1, we get that Q{u, u) is positive definite assuming that T is small. 

Having Green's formulas fl3.20p with positive definite Q('U, u) we can pro- 
ceed as in [El]. 

Let L^''\i = 1,2, be two operators of the form (12.10 and let {yo,y) = 
$j(xo, x),i = 1,2, be two maps of the form fl3.10p that transform L*^*) to L*-*) 
of the form (13.110 . i = 1,2. Let and L^^^ be two operators of the form 

(EH. 

Let vf,u{, i = 1,2, be such that L^i^u{ = 0, L^i^vf = in X^q, |y„=o = 
f,vf\y„=o = 9, u{ = ff = u{y^ = vfy^ = for ?/o = 0,2/„ > 0, i = 1,2. 
We shall denote by [fi^iolj^^^o matrices of L^*"* in {yo,y',yn) coordinates, 
i = 1, 2. As in (|3.13p we have = -g^^ , g^^^ = g^^, g~Q^ = 0. 

We assume that supp / and supp g are contained in F*^^^ x (0, T] and 
= Af ^ on r(2) X (0, T) where Af^ are the DN operators for L^'^ z = 1,2. 

Let T'f\ D^^,y}^, X^^, j = 1,2, correspond to L? , t = 1,2. It was 
proven in Lemma 2.4 in [El] that if F^^'* = F2^'* then D^q fl {yn = 0} = 
D'Sniyn = 0}. Therefore we can take Ff^ = F^ , i.e. the sets F(i),F(2) can 
be chosen the same for i = 1,2. 

o 

Denote by (YjH) the closure of C^ivj^l) in the Sobolev norm \\u\\ y(i) 

and denote by i/Q(l^-^^*]) the closure of C°° functions in Yj^^^ equal to zero on 

dYjll \ {yn = 0}. Analogously one defines H x [so,T]) and H^{T^^'> x 
[so,T]) (c.f. [El]). 

Lemma 3.1. (c.f. Lemma 3.4 in [El]) Assuming that A^^^ = A^^^ on F(2) X 
(0, T) we have 

(3.24) CiW 

for allfeHliV^^) x {.s,,T)). 

Proof: Applying the Green's formula (I3.20p for z = 1,2 and taking into 
account that A^^'' = A2^'* we get 

Q^'\u{,u{) = Q^'\ui,u{), 



21 



where Q'-*^ corresponds to L^i\i = 1,2. The inequahty fl3.24p follows from 
the ellipticity of Q^^ ,i = 1,2. □ 
Denote by Ai the domain in R"+^ bounded by the planes: T2 = {t = 
T-yn-yo = 0,0<yn<i,y' e R^-^}, = {s = yo - 2/n = 0, f < < 
T,y' e R"-^} and r4 = {yo = T,0 < yn < T,y' e R"-^}. Let Li be an 
operator of the form flS.lSp in Ai. 

Lemma 3.2. (c.f. Lemma 3.1 in [El] and Lemma 3.1 in [E3]) For any 

o o 

Vq eH^ {T2) there exists u E H^{Ai),Wo eH^ {V4),Wi E ^^2(1^4) such that 
Liu = in Ai, u\r2 = vq, u\r^ = 0, ulr^ = wq, Uyglr^ = wi. 

Proof: 

Integrating by part as in the proof of fl3.2UI) and taking into account that 
ti|r3 = we get an identity (c.f. (13.11) in [El]): 

(3.25) Q{vo,vo) = E{u,u), 

where 

/n 
(l^yol^ - di^'^vi^^k + Vi\u\'^)dy. 

Once the identity (13.251) is established, the proof of Lemma 13.21 peoceeds as 
in [El], Lemma 3.1. 

Lemma 3.3 (Density lemma), (c.f. Lemma 2.2 in [El]) For any w E 
i^o(-Rjso) there exists a sequence u-^" E HliYjso)., fn ^ Hl(T^^'^ x {sq,T)), 
such that 

llw — n-^" 111 y. ^0 when n^oo. 
Note that H^iRjso) C H^iYjJ- 

The proof of Lemma [3.31 is based on the Green's formula (I3.20p . Lemma 
I3.2l and the unique continuation theorem of Tataru (c.f. [T]) and it is identical 
to the proof of Lemma 2.2 in [El]. □ 

The main lemma used in the proof of Theorem 13.11 is the following 

Lemma 3.4. (c.f. (2.4O) in [El]) Let u{,Vi,i = 1,2, be the solutions 
of L^l^u = in X^(l,i = 1,2, with zero initial conditions and u{\y^=Q = 



f, |j;„=o = 9, where f,g belong to Hq{T''^'> x [0,T]). Suppose A^^^ = A^^ 
r^^) X (0,T). Then for any sq E [0,T] we have 

(3.26) [ ^'^^dsdy'= I ^^dsdy'. 



on 
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The proof of Lemma 13.41 uses Lemmas 13.11 and 13. 3^ and is exactly the 
same as the proof of (2.40) in [El]. We shall repeat this proof here for the 
convenience of the reader. 

Integrating by parts we obtain 

(3.27) / {uiv^ - u^vDdsdy' = 2 ! u{v^dsdy' 
Jy2o Jy2o 

uf{T,y',0)vs{T,y',0)dy'. 

Since n^(T,2/',0) = fiT,y'), v^iT,y',0) = g{T,y') we get, using dSHD, that 

{u{.,v^) jy^^ulWdsdy' is determined by the DN operator. Therefore we 
have 

(3.28) (_i,„») = (^.„|) 

for all e ^o(r(2) X (0,T)). Consider f,g e H^{T^^^ x (0,T)). Then 

supp u{ and supp ff are contained in Y^Q,i = 1,2. Take any sq G [0,T). It 

follows from the geometric assumptions that f/q^ fl {s > Sq} C -R2so' ^« 

be such that ^ = when s > sq and Wi\s=so = u{\s=so- 

Let = u{ — Wi when s > sq, u^^ = when s < sq. Assume that / and 

therefore u{ is smooth. Then e i^oH^So) ^ ^oH^i'o)- We shall prove 
that 

(3.29) {nS,vf) = {u'i],vl) 

for any g G -ffo^(r« x [0,T]). 

By Lemma 13. 31 there exists G -^^0(^250*) such that || ||-^ ^(i) 0. 



"0 

By Lemma [3?T] \\ul" — v^'^^\\ (2) for some v^'^^ G H^{Y2g^). Substituting 
f = fn (13.281) and passing to the limit when n 00 we get 

(3.30) {u^oJ,vf) = {v^^\vl). 

Note that flCTj) holds for any g G H^{T^^^ x (0,T)). Take G H^{T^^^ x 

(so,r)), i.e. G //o^(r2S)- Since 45 = ^ when s > Sq, and vf = for 
s < So, i = 1, 2 we have (c.f. fl3.28p 

(3.31) {u^^],v{) = (4?, ^f'), G ifo^(F22). 
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Compairing fl3.30l) and (13.31 1) for g = g', we obtain 

(3.32) {u'iJ,vl') = {v^^\vl'). 

Since vf G -f^o(^2so) is arbitrary we get by the Lemma [3731 that 

(3.33) = 4? on 4t 



When g G H^iV^^^ x [0,T]) we have that (supp vf) n {s > sq} C F/q ^ 

?(2) 



{■s > So} C R%so- Therefore we can replace vi^^ by u^q] in (13.301) when 



G i^oH^/o^) and this proves (K^ . Finally, substracting (IX^ from ([S2H1) 
we get dsS]). □ 
The next step of the proof of Theorem 13.11 will use the geometric optics 
solutions. Since the constructions here differ from [El], page 824, we shall 
proceed with more details. As in (2.41) in [El] we are looking for u{ in the 
form: 



(3.34) u{ = e^'^^^--) Y: jLo^i^^ r, y') + «r ^ 



TV 

where A; is a large parameter, i = 1, 2 



(3.36) *^-^]:Siv>^-0, 

ao^\y„=o = Xi{s)X2{y'), « = 1,2, 

cip\p > 1, satisfy nonhomogeneous equations of the form (13.351) that we will 
not write here and u^^^^^ is the same as in (2.41) in [El] (c.f. [El], page 
824). Here Xi(s) G C^iR^), Xi(s) = 1 for |s - So| < (5, Xi(s) = for 
\s — so\ > 26, S is small, X2{y') £ C^(r'^-'^)) is arbitrary. 

Let f3j^\yn,a) be the solution of the system of differential equations 

(3.36) -^ = 2(7°J(/5«,2/„), /3f\0,a)=aj, 1 < j < n - 1, t = 1,2. 
dyn 

Let a^*) = {a^^\yn,y')} be the inverse to = {l3j^\yn,a)}- We have 

(3.37) 9°i"(y.^') _g^g(y)||!^o. «fU.„ = a. l<j<n-l. 

k=l 
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Therefore aQ{s,T,y') = Xiis)X2{ct^^\'^—^—^,y')) is the solution of fl3.35p . 

aQ''|j;„=o = Xi{^)X2{y')- Substituting the geometric optics solutions (13.341) 
in fl3.26p . integrating by parts and taking the limit when A; ^ oo we obtain 
(c.f. (2.42) in [El]): 
(3.38) 

/ X2{a^'\^^^,y')vfiso,0,y')dy' = [ X2{a^'\^^.y')vl{sM)dy'- 

Note that r = on Y^^ ,i = 1,2. Changing T to T - r', < r' < T we get 
(I3.38P for any < r < T. Consider the following change of coordinates 

(3.39) s = s, T = T, = a«(^-l|^,ty), ^ = 1,2. 
The inverse change of variables has the form: 

(3.40) S = S, T = T, y' = p(^)(l^l^—^._y'), z = 1, 2. 

Note that y' = f^^^\yn,y') is the endpoint of the curve (13.361) starting at 
y' G r*-^^ when yn = and y' = a^'\yn, y')- 

Let S = {(s, r) : s > 0, T > 0, s + T < T}. Denote by /9(^)(S x f^^^) the 
image of S x f^^^ under the map ([330]), z = 1,2. Note that (3^^{T. x f^^^) 
is contained in x[q. Therefore Ril Qi fl (3^^\Tj x T^^^) is contained in 
RiQ,i = 1,2. Here Qi is the rectangle {{s,T,y') : r = 0, s G [0,T],y' G T^^^}. 
Denote by the image of R^q under the map (I3.39p . Finally, denote by S^*) 
the projection of R^q on the plane y^ = 0. Note that B^'^^ C T^^^ x [0, -j], z = 
1, 2. We shall assume that B^^^ D F x [0,^]- This assumption always can be 
satisfied when T is small enough. 

Make the change of variables (13.401) in (13.381) . We get 

(3.41) / X2{y')vfis, r, y'))Myn, y')dy' 

= [ X2{y')vl{s,r,P^'\ ^ ~l~ \ y'))Uyn,y')dy'. 

yn = "^--f--^5 Ji (^~2~^ ^i)') is the Jacobian of the map (I3.40p . Since X2{y') is 

any C^(F''^'*) function we get that for any y' G F*'^'* 

(3.42) 

vf{s, r, y'))Myn, y') = vl{s, r, y))J2{yn, y')- 
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Note that (Km holds for {s, r, y') e J: x r'-^\ 

Let Xi(s) be the same as before, and Xsiu') ^ C^(r'^^)) be arbitrary. 
Construct v^f^ as geometric optics solution (13.341) with g = Xi(s)x3(y')- 

Take s = so and A; — > oo. We get 

(3.43) vL = xi{so)X3{cy'^'\so,T,y')), 

where v^^ = limjt^oo fc- Substituting v^f^ in (I3.38P and taking the limit 
when /c oo we obtain 

X2{a^^\so,r,y')x3{a^^\so,T,y')dy' = / X2{a^^\so,T,y')x3{a^^\so,T,y')dy' . 

Make the change of variables (13.401) . Since X2, Xs ^-^e arbitrary we get, as in 
([322D, that Ji{y) = J2{y). Therefore 

(3.44) vl{s,T,P^'\ ^~l~\ y'))=vl{s,T,(3^'\ ^~l~\ y')l 

where {s,T,y') G S x r^^\ 

Let wf{s, r, y') = (s, r, /^W), i = 1, 2. Then (s, r, y') = u;f (s, r, y'), V(s, r, y') G 
S X T^'\ 

Our strategy to complete the proof of Theorem 13.11 will be the following: 
Making the changes of variables (13.401) in L^ivf = we get L^iwf = 
0,2 = 1,2. Using that wf = wl for all g G HliV^^^ x (0,T)) and using the 
density lemma 13.41 we shall prove that the coefficients of L^^'* and are 
equal. Since the density property holds for r fixed we have to take care of 
terms in Lf^ that contain derivatives in r. □ 
Note that integrating by parts as in (I3.27P we get 

{u{v^ — u^Vg)dsdy' 

Y20 



u^vldsdy + I u^{T,y',0)v^{T,y',0)dy' 
Jr(i) 



'1^20 

Therefore as in (13.281) we conclude that 

(3.45) iu{,vl) = iuivi). 

Using ( ]3.45p instead of ( 13.28P we get an equality of the form ( 13.26^ with the 
roles of u-^ and reversed: 
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(3.46) / ui^dy'ds = [ ui^dy'ds. 
From fl3.46p we get, analogously to (13.441) . that 

(3.47) = onSxr . 

We used here again that Ji{yn,y') = J2{yn,y') in S x r^^\ Differentiating 
wf{s,T,y') = vf{s,T,(3^) in s and y' we get 

^'•''^ — d; — = — d; — + ^ — — 



(3.49) 



dwf{s,T, y') 



n-l 

E 

fc=l 



dyk dyj 



where (5^^ = (3^^{y n,y'), yn = 
It follows from (Km that 

(3.50) 



T~s- 



dv^is, T, /?»_) ^ g dafi^, /3«) dwfis, r, y') 



dyk 



k=l 



dyk 



where 



dyk 



is the inverse matrix to 



dyj 



Substituting (l33nll into (EJHD, using ([331, we get 

(3.51) 

^ a«f)(j/„,/jW) (,) a^^f(,,r,yO ^ ^ daf\y^,P'^')) ,,^ dwl{s,r.,y') 
f^, dy^ dm f^, dyk dy, 



2:-f-^,r = o,(5,y')er^'^x[o,r]. 



where yr, 

Since {vl{s,T,y'),g G C^iV'^^'^ x (0,T)])} are dense in (R^^J) (c.f. 

o 

LemmalXSD, we get that {wf{s,T,y')} are dense in (Rio)^ where h!^q is 
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the imaere of R^^d C i?[o^ under the map fl3.38p . Therefore we get (c.f. the 
end of $ 2 in [E3]) that 

on ^[q\ Here r = 0,y„ = Note that S^^^ is the projection of ^[q'* on 
the plane Hq = 0. Therefore (13.521) holds on 5*^^^ since a*^*^ and do not 
depend on yo. We have on S x f^^^ (c.f. (Km . Km ): 

(3.53) c/f( ^~^'\ /?«(^-l|^, y')) =y^,l<j<n-l, z = l,2. 
Differentiating ( 13.53^ in s we get: 

(3.54) ««(^:i|^,/3«(^^^|^,2;0) 

n— 1 

+ E (^^f^, y )) = 0, . = 1,2. 

fc=i 

Combining (I3.54p and (I3.52p we get 
(3.55) 

a'lhyn,P^'\yn,y')) = a'^hyn,P^'\yn,y')), l<J<n-l, {y^,y') G 

Consider the equations L^ivf = in X^f^. It has the following form in 
{s,T,y') coordinates: 



where g'^^^ = g^. Note that g^^^ , i.e. the coefficient of , is zero. 
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Making the change of variables fl3.39l) we get equations of the form 

(3.5r) Vf^^.,M{l.,S^,l.j^ 



«(y„,/3«)«;f(s,r,y') = 0, (s,r,^') G S x V^'\ 



where 
(3.58) 

ivn, y)=l^^ gfo ivn. 

1 < j, A; < n - 1. We used in fl337D that (c.f. (EST])) 

dyp dr 



(3.5.) a-^i:.^'%../^-') ^°-"'"f'""'" -^-°. 

p=i 



daf{y,^,P^%y^,y')) . , , ^ . ^ , 

^io = ' since =0, l<J<n-l. 

Since wf^s, r, = ti^f (s, r, y') in S x r*'"^\ we have in B^^^: 



(3.60) (L«-Lf))^?=$:j-A(j, 



(^10 - ^2o) 



dyk 

+{vi'\yn.P^''')-Vi^\y^,(3^^^))wl = Q. 
We took into account that Ji (?/„,?/') = J2{yn,y') holds on F*'"'^^ x [0, -j] and 

«5^(?/n,/^('Hyn,y')) = «5?(?/„,/?(')(l/n,y')), 1 < J < ^ " 1, holds OU 

o 

Since {wf,g e Co"^(r(i) x (0,r))} are dense in (R^^J), we get, as in 
[E3] (see the end of section 2 in [E3]), that 

(3.61) ~gto = ~9i^ Vl'\yr„f3^''>) = vl'\y^,p'^'^) in ^ff. 
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Noting that the coefficients in fl3.6ip do not depend on uq and B^^^ is the 
projection of R.'^J on uq = we have that (13.611] holds in B'^^\ Therefore we 
proved that L^-^^ = Lf'^ in B^^K Now we shall prove that also L^^) = L^^) m 
B^^\ where L^*-* is the opeators L^^^ (see (13.81) ) in {s,T,y') coordinates. 
Operators L*^*-* have the following form (c.f. (13.571) ): 

(3.62) L« ^ (|^--(.„,/3«(.„,^'))^ + ^J^^-^y-^''^^s 



where gj'^ has the form (I3.58P with gf^ replaced by gf^iun, P'^^Kun, y'))- Since 
9Z = {at^'Y^gT we get that 

(3.63) ~gl\yn,y) = {9p~iyn,P^'^)r'~gloiyn,y)- 

We used in ( Km that g~'^ = for 1 < j < n - 1, z = 1, 2, and that (Km 
implies 

ay. ar ' 

j=l 



since = p^. 

Therefore to prove that L^^^ = L^"^^ it remains to prove that 

(3.64) ^+-(y„,/9«)=^+--(y„,/9(2)). 

Making the change of coordinates (I3.39P in (I3.14p we get 
(3.65) 

K«(. /?«(. v'))--S-J-^^(ja^'^]-Ta^''^^ 
V, (y„,/3 (,.,,))- yj.g^, J ^2^^,.o 5^^. , 

where ?/„ = ?/„, iW(?/n,?/') = v4W(?/„, /?W(?/„, ^')). ^io"' 1 < i, ^ < - 1, are 
the same as in p.57p, ^^q" = ^lo = -«!?(?/« , 1 < i < - 1, ^IT = -1- 
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Taking into account that qIq = g^Q, Ji = J2, and 



Uj Uk iJj yk Uj yj ^ yk ^ yk yk ' yj 

we can rewrite 



as homogeneous second order elhptic equation for A^^^ (3^^^)—A^'^^ {yn, /S*-^-'), 
where y') = ln(a+'-(y))i = \n{j^{det[gt{y)]];,U)--^ (^-f- O). 

Since A^^^ and A^'^^ have the same Cauchy data when ?/„ = (see Remark 2.2 
in [El]) we get, by the unique continuation theorem for the eUiptic equations, 
that A(i)(y„, /?(!)) = A(2)(|/„,/3(2)) in 5(1). Since 



9i 



and 

we get fl3.64p . Therefore L^^^ = iP"^ in B^^K Note that by the assumption 
i?WDrx[0,f]. □ 

Theorem 13.11 concludes the local step of the proof of the main Theorem 
12. 3[ The global step of the proof is similar to the proof in [E2] : 

Consider the initial-boundary value problems for L^^^Ui = in domains 
fiW = fio \ Uj"^i^5'\ i = l,2. Let Ai C be the image of v?"^ o a^^{T x 
[0, -j]), where ctj is the map (13.391) and $4(2:0, x) = (a;o + ai{x) , (fi{x)) is the 
map (13.101) . Denote by $3 the map $3 — $-|^ ^ o G;^^-^ o ^2 o where P2 has the 
form (13.401) . Note that $3 is a diffeomorphism of the form (13.101) . $3 = / on 
(A2 n dQo) X (—00, +00) and 

$3oL(2)=L(^) on Ai. 

Note that any map $ of the form (13.101) can be represented as a composition 
^ = ai o tpi = {p2 ° 0,2, where ipi are the diffeomorphisms of A2 onto Ai and 
maps tti have the form yo = xq + aj(x), y = x, ai{x) G C°°, ai{x) = on 

dno. 

It follows from [Hi], Chapter 8, that there exists an extension $3 of the 
map $3 such that $3|ef7ox(-oo,oo) = I, $3 has a form ( 13.101 ). i.e. $3 = 03 o 

ips, LPs is a diffeomorphism of fi*"^"* onto Q^^^ ips{Q'''^^). Denote L^^^ = as o 
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^3oL^^\ Then L^^) is a differential operator of the form fl2.ip on Ai C 
and L(3) = L(i) on Ai. 
The proof of the following lemma is the same as in [El], Lemma 3.3 (c.f. 
[KKLl], Lemma 9): 

Lemma 3.5. Let A[ C Ai be such that Oi \ A^ has a smooth boundary, 71 = 
dflo n dA[ is connected and L^^^ = L^^^ on A[. Let 72 = dA[ \ 71. Suppose 
A^^) = A*^^) on dilo X {—00, +00), where A^'^^ are DN operators corresponding 
to LW, respectively, i = 1,2,3. Note that A^^^ = A*^^) on d^o x (—00, +0x3). 
Then the DN operators A^^^ , A^^'* corresponding to the operators L'^^^ , L*^^-' in 
the smaller domains Vt'^'^'> \ A-^^Vt'^'^^ \ A^ are equal on {{BVLq \ 7^^) U 72) x 
(-00, +00). 

Therefore Theorem 13.11 and Lemma 13.51 reduce the inverse problem for 
/^(i)^/^(2) in ^(1) X (— c)0,+oo), QS"^^ X (—00, +00) to the inverse problem for 
/,(3) in smaller domains (fi(^)\Ai) x (-00, +00), a'J x (-00, +00). 

Continuing this process as in [E2] we can prove the main Theorem 12.31 
Note that it is enough to have A*^^-* = A*^^) on i7 x (0,To), where Tq is large 
enough, to prove Theorem 12.31 
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